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The family of fifth order nonlinear evolution equations is studied. Some traveling 
wave elliptic solutions are found. The classification of these exact solutions is given. 
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1 Introduction 
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Recently in Refs. [HE] a method was introduced for classification of meromor- 
\ phic exact solutions of nonlinear ordinary differential equations. This method 

allows us to classify and construct meromorphic exact solutions for a wide 
class of autonomous nonlinear ordinary differential equations in the explicit 
form. The main idea of the approach from [HE] is to compare the Laurent 
series corresponding to solutions of ordinary differential equations with the 
Laurent series for the general form of possible meromorphic exact solution. 
The main advantage of the method from works [HE] is that we can construct 
and classify more general forms of exact solutions in comparison to existing 
methods for finding exact solutions P-fl5]. Using this approach in the work p], 
meromorphic exact solutions of second-order differential equation were clas- 
sified. The authors of [3] presented the classification of elliptic solutions of 
a third-order differential equation obtained with the help of the above men- 
tion method. This method was used for constructing and classifying traveling 
wave solutions of the Kawahara equation [2] and the generalized Bretherton 
equation [TH] as well. 
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In this work we study traveling wave solutions of the equation that takes the 
form 

u t + au m u x + fiu xx + j3 u xxx + vu xxxx + 5u xxxxx = 0. (1.1) 

Eq. (11.11) at m = 1 was obtained in [T7] for the description of nonlinear waves 
in a viscoelastic tube. Nonlinear wave processes described by (II. ip at m = 1 
were studied numerically in [IB]. The meromorphic solutions of Eq. (II. ip at 
a = 7 = were found and classified in [2J. Elliptic solutions of Eq.f ll.ip at 
/i = /3 = ^ = 0in the cases m = 1 and m = 3 in terms of the Jacobi elliptic 
function were found in [7]. In [7] elliptic solutions of Eq. fll.ip in the case of 
\i — v — and m = 2 were obtained as well. The simple periodic solutions of 
Eq. (Il.ip at any m were found in work [19] . However elliptic solutions of Eq. 
(II. ip at v 7^ and /i ^ were not considered previously. 

Using the traveling wave u(x,t) = y(z), z = x — C t in (11.11) and integrating 
the result we obtain 

Ci - C y + —y m+1 + fiy z + (3y zz + uy zzz + 5y zzzz = (1.2) 
m 

Here C\ is an integration constant. 

The aim of this work is to construct and classify elliptic solutions of Eq. (ll.2p 
at m — 1, 2, 4. To this aim we use the approach suggested in p~H3]. 

This paper is organized as follows. In section 2 we give a brief description of 
the method from [IHS] and present elliptic solutions of Eq. (ll.2p at m = 1. 
In sections 3 and 4 we construct elliptic solutions of Eq. (ll.2p in the cases of 
m = 2 and m = 4 respectively. In the last section we discuss our results. 



2 Meromorphic solutions of the equation studied at m — 1. 

Let us briefly describe the approach from works [IHH] for finding elliptic solu- 
tions of nonlinear ordinary differential equations. 

We assume that solutions of Eq. (11.21) can be presented in the form of Laurent 
series in a neighborhood of a movable pole z = z 

oo 

v(*) ro = £a?(*-*>)*- , '> P,ieN (2.1) 

fc=0 

In the case of m = 1 solutions of Eq. (ll.2p admit one expansion (Z = 1) in the 
form (12.11) . At m = 2 and m = 4 solutions of Eq. (11.21) admit two (/ = 2) and 
four (/ = 4) expansions in the form (12. ip respectively. 
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Eq.f ll.2p is autonomous and thus without loss of generality we can set zo = 0. 
Using the residue theorem for an elliptic function we obtain the necessary 
condition for existence of the elliptic solutions. This condition is the follow- 
ing: the sum of coefficients of z~ l in series ( 12.1 j) is zero. From the necessary 
condition follows that while Eq.f ll.2j) possesses one expansion of type (12. ip 
the coefficient at z~ x has to be zero. In the general case we can present this 
necessary condition in the form 

X)4-i = ' « = 1,2,3,... (2.2) 
i 

From here and below in the case of one expansion of type (12 .1|) we omit the 
index I in coefficients of the series. 

We use the method from p~H3] for constructing the elliptic solutions for Eq. fjl.2j) . 
The algorithm of this method is the following: 

(1) Construct the formal Laurent expansion of type (12. ip for a solution of 
Eq.dOD; 

(2) Check the necessary condition (12. 2 j) for existence of the elliptic solutions; 

(3) Take the general form of the possible elliptic solution for Eq.( ll.2p and 
find the Laurent expansion for the possible solution; 

(4) Compare the formal Laurent expansion of solution for Eq. (11.21) that was 
found in the first step with the formal Laurent expansion of the possible 
elliptic solution that was found in the third step; 

(5) Solve the system of the algebraic equations obtained in the fourth step 
and find the parameters of Eq.f ll.2j) and parameter of the possible elliptic 
solution. 

We can see that the second step of the algorithm coincides with the first two 
steps of the Painleve test. The Laurent expansion for possible elliptic solution 
can be find by using textbook [20] or by using symbolic computation software. 

Let us construct elliptic solutions of Eq. fll.2p at m = 1. Without loss of gen- 
erality we assume that a = 6, 5 = —1 in Eq. fll.2p . Thus from Eq. fjl.2p we 
have 

Ci - C y + 3y 2 + fi y z + (3 y zz + v y zzz - y zzzz = (2.3) 

Eq.f l2.3p admits one Laurent expansion (I = 1) of type f|2.ip in a neighborhood 
of moveable fourth order (p = 4) pole. The Fuchs indices corresponding to 
expansion in the Laurent series are the following 

ji = -l, .72 = 12, j 3>4 =i(ll± 4 yi59) (2.4) 
We see that the Fuchs index j 2 has a positive integer value. Thus the expansion 



3 



for the solution of Eg. (12. 3 j) can exist if a\i is an arbitrary constant. 

The expansion for the solution of Eg. ( 12. 3 p in the Laurent series is the following 



280 280z/ /21z/ 2 7f3\ 20 / 59i/ 3 20/3// /A 140 
~ + Tlz 1 ~~ [~m + YJ ilk 2 + ^17303 + 1287 + 9 J 23~I + (2.5) 

+ . . . + a i2 z 8 + ... 



Series (12.51) corresponding to solution of Eg. (l2.3p exists in the case 
784354/i 2 /3 2 z/ 2 + 137745805669/^ /3 V 449029786141 1/x z/ 9 



3671430867 549479397531636 297960291083728988 
35385728065 1/jV 2781278216285z/ 10 /3 206626426493z/ 6 /3 3 + 



4029515581898664 297960291083728988 3296876385189816 
1673496887869z/ 8 5 2 310/3V 2 194170261339/i v 7 
+ 44324671400885304 + 65033397 + 1704795053880204 + 
, 163837/i 3 /3 i/ 32794140175/z 2 /?z/ 4 80582613855300z/ 12 



1498986567 99905345005752 99146286858110820757 
(5313276^ + 7643363/x 2 ) (10286502336a + 25164119899k 3 



506657459646 674361078788826 
v 



85961304 

jw (42262166548/3 3 + 40586709267/i 2N ) = 

1348722157577652 V ' > 



(4961 v/3 + 19844 // + 10737 z/ 3 ) (C 2 - 12(7 X ) + 



(2.6) 



The last eguality is the compatibility condition for existence of the Laurent 
series (12. 5p . Series (12. 5 p does not exist if relation (12. 6p is not satisfied. 

In accordance with the classification of meromorphic solutions of autonomous 
ordinary differential eguations presented in [H|2] there is one type of possi- 
ble elliptic solution of Eg. (I2.3j) . This are elliptic solutions corresponding to 
Laurent series (12. 5p . 

Let us construct elliptic solutions of Eg. (12. 3 j) . Taking into account the neces- 
sary condition for elliptic solutions to exist ( 12. 2 j) from series ( 12. 5 j) we obtain 

_ (531z/ 2 + 2420/?);/ , 2 ^ 

^ ~ 17303 

Using the approach from [HE] we find that the general form of the possible 
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elliptic solution of Eq. fl2.3p corresponding to series (12. 5p has the form 
20 (21u 2 ^ 7/5\ UOu 

y = "is [ui + T) p{z > 9l > 92) ~ TT P ( ^^ 2)+ ^ 

+280 p 2 (z,g u g 2 )-^ + h 



Comparing expansion (12.51) with the Laurent series for (12.81) we obtain 
C 1 /4845z/ 7/3 2 643£z/ 2 \ 

tin = 1 1 , 

6 169 V 29282 9 726 / 

(121/3 + 48^ 2 ) (2057,5 + 504z/ 2 ) / 2 .9) 
92 ~ 207843636 ' 

_ v 2 (769365 iA + 6912246 (3u 2 + 19399325 (3 2 ) _ 41/3 3 
93 " 934108684080 2372760 

Taking into account (12.71) from condition (I2.6P we have 
3u 2 



(2.10) 



C x = ( 17324377875 u 6 + 187507377750 u 4 p+ 

171424512006748 V 

+709528086443 /?V + 1078179110844 (3 3 ) + - 1860/3 
H ' ) 12 199927 

Also using (12. 7p and (I2.10p we obtain the value of a\i 

/3 5 (51717941/3 + 174854187// 2 ) 27z/ 10 (322959854223z/ 2 + 5492718054164/3) 



1272898705720561920 38819181898104129541859840 
i/ 4 /3 2 (14947544754846/3z/ 2 + 4344006788 lOOi^ 4 + 26834189262353/3 2 ) 
167038348444816621893120 

(2.11) 

Substituting h , g 2 and (12 .9p into formulae (I2.8P we obtain the elliptic 
solution of Eq.Q. 

Let us note that at (3 = — and /3 = — 3 j 2 y elliptic solution (12 .8p degenerates 
to simple periodic solutions. 



3 Meromorphic solutions of the equation studied at m=2. 



Let us consider the traveling wave solutions of the equation at m — 2. From 
Eq. (11.21) we have 

C 1 -C y+ -y 3 + ny z + (3y zz + v y zzz + 5 y zzzz = (3.1) 
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Without loss of generality we assume that 5 = — 1 and a = 360 in Eq. (13.1j) . 
Thus from Eq. (13. 1 1) we obtain 

C x - C y + 120 y 3 + \i y z + (3 y zz + v y zzz - y zzzz = (3.2) 



Eq. (13.21) admits two different Laurent expansions (/ = 2) in a neighborhood 
of a moveable second order (p = 2) pole. The Fuchs indices corresponding to 
expansions of the solution in the Laurent series are the following 

j'i = -l, j 2 = 8, j 3 ,4 = ^(7±2v / 7l) (3.3) 



We see that the Fuchs index j 2 has a positive integer value. So expansions for 
the solution of Eq. (l3.2p can exist if a§ is an arbitrary constant. 

Expansions for solutions of Eq. (13. 2p in a neighborhood of a moveable second 
order pole are the following 

y M> = ± l ± ^ T 23^ + 98^ 1 / M Jhi M , + ... + Mi . + .. . 
y z 2 Uz 5880 4 V 45 1715 126 8 



(3.4) 



Series (13 Ah corresponding to solutions of Eq. (I3.2D exist in the case 
, z/ 2 /3 3 ± /V + 3u 2 d ± 2323z/ 3 // ^ 2C ^^ ^ 31C u A , 



"2100 900 28 617400 315 27440 ~ 

577^ 8 C Q u 2 (3 , 1289z/ 4 /3 2 , 31z/V , 257^ 6 /3 , , Q ,x 

zp — - — H± ^_ ± C— ± c_± (3.5) 



"9882516 630 " 1234800 8400 540225" 
, 4381 u 5 fi ± nP 2 v _ 



"4321800 630 



In accordance with the classification of meromorphic solutions of autonomous 
ordinary differential equations presented in |1]|2] there are different types of 
possible elliptic solutions of Eq. (13. 2p . The first type is elliptic solutions corre- 
sponding to one of the Laurent series ( 13. 4ft . The second type is elliptic solutions 
corresponding to both of the Laurent series ( 13.41) . 

The necessary condition for existence of elliptic solutions ( 12. 2 j) in the case of 
one of the Laurent series (13 ,4p gives us v = 0. In the case of elliptic solutions 
corresponding to both of the Laurent series (13.41) the necessary condition (12. 2p 
is automatically satisfied. 

Let us construct elliptic solutions of Eq. (I3.2j) corresponding to both of the 
series (13 ,4p . In this case we force that the compatibility conditions (13. 5p are 
satisfied simultaneously. We see that this is correct at C\ = or at v = 0. 
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In accordance with the method presented in [TJ[2] the possible elliptic solution 
of Eq. (13. 2p corresponding to both of the series (13. 4 p has the form 



y = - 



p'{z,g 1 ,g 2 ) + B 
p(z, 91,92) - A 



v P {z,g u g 2 ) + B ,„ R x 

™~ 7 v - + 2p{z,gi,g 2 ) + h (3.6) 

28 p(z,g 1 ,g 2 ) - A 



Here A, B and h are parameters that will be found later. Also we denote 4- 
by'. 

First we consider the case of G\ — 0. Comparing expansion (13.41) with the 
Laurent series for H 3 . 6 1) we obtain 



„ v 3 - 196u n n 38416 a 2 + z/ 6 - 532 u 3 a 
h = A= — — r , B = 0, g 2 = 



#3 



1680// 

196 n - u 3 ) (38416 /i 2 - 567 u 3 fx + u 6 
296352000// 3 



141120z/ 2 



P 



39u 3 - 1372// 
196^ 



(3-7) 



From the compatibility conditions (I3.5P we have 



134456 /i 2 - 4067 z/ 3 a + 11 z/ 6 
°° ~ 6860^ (3 ' 8) 



Solving the algebraic system of equations for parameters ho, A, B, g 2 , g 3 we 
obtain values for constant as in the form 



1 

23897825280000z/ 



4 1,2) = ooonToo C oonnnn. 4 f =•= 1475789056 ^ 4 T 8605184 i/VT 



(3.9) 

^263424 z/V 2 =F 224 z^V ± f 1 



Expressions (13. 9 j) are necessary conditions for existence of the Laurent series 
(13. 4p . These conditions show us that the elliptic solution of Eq. f]3. 21) contains 
only one arbitrary constant. In this case we can add the arbitrary constant zq 
to the variable z. 

Taking into account formulae (I3.7P we have an elliptic solution of Eq. (13.21) in 
the form 



p'(z : gi,g 2 ) 



vp (z, 9i, &) 



+ 2p(z, 91,92) + ho 
(3.10) 



2S(p(z, 9l ,g 2 )-^^ 
where g 2 , #3 are defined by (13.71) . 

The elliptic solution (13. 6 j) degenerates to simple periodic solutions at the fol 
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lowing values of fi 

v 3 \lv 3 v 3 v 3 

^= 476' ^= 1372' ^= 686' ^=56 (3 " U) 

3 

In the case of /i = ^ we obtain the following simple periodic solution from 
solution (13. 6p 

= ^(7-14cos^{^} + v^T9sin{^}) 
V ~ 13328 cos 2 { } sin 2 { } 

In this case parameters A, ho, 8, Co can be written as 

h = A = , 8 = , C = (3.13) 

2856 833 1586032 



At fj, — solution (13. 6p degenerates to 

y= V - r 4 (3.14) 
y 49(e^-l) 2 

Parameters A, ho, 8, Co are defined by the following relations 

V 2 a 21/2 n 15z/4 / Q1 ^ 

/i0 = A= -TT76' ^ = -- Co = -5488 (3 - 15) 



In the case of ji = ^ from (I3.6P we have 

y = ^ + , 2e " ~ 1} (3.16) 

^ 392(eif-l) 2 1 ; 

In this case parameters A, ho, 8, Co can be written as 

y2 a 41z/2 r, 151/4 Cx^\ 

h = A = , 8 = , C = (3-17) 

2352 H 196 19208 



3 

At fi = |g solution (I3.6p degenerates to 



56 

V 2 



7-sin{fe}) 



1 568 cos 2 {4|p} 
Parameters A, ho, 8, Co are defined by the following relations 



(3.18) 



= = -™*, C =-i^ (3.19) 

672' ' 392 ' 5488 
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In the case of v = 0, the elliptic solution of type (13 .6p degenerates to a simple 
periodic as well. 



Let us construct elliptic solutions of equation (13 .2p that correspond to one of 
the series (13 .4p . We see from (13 .4p that the necessary condition for existence 
of an elliptic solution is v = 0. In this case the possible elliptic solution of Eq. 
(13. 2 p corresponding to the first series from (13 .4p is the following 

V = p(z, 92,93) + h (3.20) 

Comparing the Laurent expansion for (I3.20p with the first expansion (13. 4 p at 
v = we have 



, P _ IOCq - (3 2 5C /?-450C 1 -/3 3 , , 

fin = , U = 0, flo = , 03 = [o.ALj 

° 60' P y 180 ' y 5400 



The coefficient a§ is defined by the following relation 

(10C - /3 2 ) 2 



a 8 



38880000 



(3.22) 



Substituting (72 , ^3 and ho from (I3.2ip into (I3.20p we obtain the elliptic solution 
of Eq.(GL2D at fi = v = 

/ IOCq - /? 2 5^0/3-450^-/3^ _ ^ 
^ P \ ' 180 ' 5400 / 60 1 ' ' 

This solution was previously obtained in work [2]. 

Thus we see that equation (I3.2p has elliptic solutions corresponding to one of 
the series (13 .4p only at v = fx — 0. Elliptic solutions of Eq. (l3.2p at v = \i — 
(3 = can be obtained form the formulae (I3.23p . 

Note that Eq. (l3.2p possesses the symmetry 

y(z,-C 1 ) = -y(z,C 1 ) (3.24) 

Thus the elliptic solution corresponding to the second series from (13.41) can be 
obtained from solution (13.231) . 
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4 Meromorphic solutions of the equation studied at m=4. 



Let us construct elliptic solutions of Eq.f ll.2j) at m — 4. In this case from 
Eq. fll.ip we obtain 



C 1 -C y + ^y 5 + Liy z + Pyzz + v y zzz + 8 Vzzzz = (4.1) 
5 



Without loss of generality we assume that 5 = — 1 and a = 120 in Eq.f l4.ip . 
Thus from Eq.f l4.ip we have 

C\ - C y + 24 y 5 + fx y z + (3 y zz + v y zzz - y zzzz = (4.2) 



Eq. (14. 2 ft possesses four different expansions in the Laurent series in a neigh- 
borhood of the first order moveable pole. The Fuchs indices corresponding to 
these expansions are the following 

ii = -1, J2 = Q, j 3 ,4 = ^(5±zv / 39) (4.3) 

We see that the coefficient a§ in the Laurent series has to be arbitrary. Oth- 
erwise the Laurent expansions do not exist. 

Expansions of solution of Eq.f l4.2j) in a neighborhood of the moveable pole are 
the following 

2/ (1 < 2) = ±\ ± |j T (3^ 2 + 10/3) Z + ... + 4 1,2) z 5 + • • • (4-4) 

y {3A) = ±- z ± | + ^ (3^ 2 + 10/3) ,+ ...+ 4 X ' 2) , 5 + . . . (4.5) 

The compatibility conditions for series ( 14. 4p to exist are the following 
^ 189i/ 6 z//3 (10z//3 + 35/i + 9 u 3 ) fi (99z/ 3 + 100 fji) 



100000 750 " 3000 (A6) 
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In the case of series (14. 5 p the compatibility conditions have the form 
, 189ii/ 6 zz//3 (10^/3 + 35/^ + 9 z/ 3 ) i> (99z/ 3 + 100 



100000 750 3000 ^ ^ 

(— 5zCi ± Co^) 

id = 
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Using the necessary condition ( 12. 2 j) for elliptic solutions to exist we see that 
Eq. (I4.2j) may admit three types of elliptic solutions. There are elliptic solutions 
corresponding to both of the series (14.41) or to both of the series ( I4.5p . Also 
Eq.( )4.2p may admit elliptic solutions corresponding to series f !4.4j) together 
with series (14. 5p . 

Let us construct elliptic solutions of Eq. (l4.2p corresponding to expansions 
( I4.4p . The general form of possible elliptic solution of Eq.( l4.2p is the following 

w 

p'{z,g 2 ,g 3 ) + B 
2{p{z,92,g3) - A) 

Compatibility conditions H4.6I) are the same at C\ = or at v = 0. In the 
case of C\ = 0, the algebraic system of equations for parameters (3, v, Co 
is inconsistent. Thus we have to consider the case of v = 0. Comparing the 
Laurent expansion of (14. 8 p with series (I4.4p we find 



H = h = B = d = 0, A = -—, g 2 



93 



60' 120 
P(13[3 2 - 150C ) 



(4.9) 



108000 

Coefficients a^' 2 ^ are defined by the relations 



fi,2) = (3(l8C -/3 2 ) (41Q) 
6 302400 



Substituting A,B,h from (14 .9p into formulae (14 .8p we obtain an elliptic so- 
lution of Eq.flOU at n = v = C x = 

w = - f P(*'a»f ) . (4.11) 

where g 2 , gs are defined by (14.91) . Solution (14.111) was obtained in [2]. 

In the case of series (14.51) . the compatibility conditions (14. 7p are the same at 
C\ = or at v = as well. In the case of C\ — algebraic system of equations 
for parameters /U, f3, is, C is inconsistent again. For obtaining elliptic solutions 
corresponding to series (14. 5p we take into account the following symmetries of 
Eq.flOJ 

y(z,-C 1 ) = -y(z,C 1 ) i y{z,iC x ) = -iy{z,C x ) (4.12) 



Using (jOl) and (14.111) we have 



ip (z,92,g3) 



w = — (4.13) 
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where g2, <?3 are defined by (14. 9[) again. Coefficients ajf' 4 '' are given by the 
following relations 

(3,4) = Waco-?) (414) 

6 ^ 302400 

In the case of elliptic solutions with poles of four different type, the algebraic 
system for the parameters is inconsistent. 

We can consider the case of /i = /3 = v = 0. However while this condition is 
satisfied elliptic solutions of type ( 14. 8 p degenerate to simple periodic solutions. 



5 Conclusion 



In this paper we have studied elliptic traveling wave solutions for the family 
of fifth order nonlinear evolution equations. We have presented a classification 
of elliptic solutions of this family in the cases of m = 1, m — 2 and m — 4. 
The explicit form of the elliptic solutions is given as well. In the case of m = 1 
we have obtained a new elliptic solution with one pole in the parallelogram of 
periods. At m = 2 a new elliptic solution with two poles in the parallelogram 
of periods is presented. We have shown that elliptic solution with one pole in 
the parallelogram of periods exists only in the case v = \i = 0. In the case of 
m = 4 we have given elliptic solutions with two poles in the parallelogram of 
periods. As well, for elliptic solutions with four poles in the parallelogram of 
periods the corresponding algebraic system is inconsistent. 
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